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ABSTRACT 


A GPSPAC/Landsat-D Interface (GLI) Ground Support System was 
built to validate the performance and to calibrate the accuracy 
of the experimental navigation package, GPSPAC, flown on the 
Landsat-4 and 5 spacecraft. Although the GLI system operated 
successfully to give the orbit information needed to validate the 
GPSPAC, it also detected two anomalies: one is characteristic of 

the GLI system and the other is characteristic of the pre- 
operational phase of GPS. Several methods were applied to 
resolve or reduce the anomalies. This paper presents a 
description of the problems, the methods applied to resolve or 
reduce them, and the results. 
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2.0 GLOBAL POSITIONING SYSTEM BACKGROUND 


The Global Positioning System (GPS) is an advanced satellite- 
based navigation system, being deployed by the Department of 
Defense, that will provide extremely accurate position, velocity, 
and time information to a variety of users 24 hours a day, 7 days 
a week. Both the Landsat-4 and 5 spacecraft carried an 
experimental navigation package, the Global Positioning System 
Package (GPSPAC) , to assess the performance and the accuracy of 
the onboard use of GPS data. 


The GPS configuration consists of a Master Control Station (MCS) 
and a constellation of Navigation Development Satellites (NDSs) . 
In its operational configuration, the NDS constellation will 
consist of 18 Space Vehicles (SVs) in six nearly circular orbits 
of 12-hour periods (20,200 km altitude) each inclined 55 degrees 
to the equator. However, when Landsat-5 was launched in March 
1984, the NDS constellation consisted of five operating SVs in 
two orbit planes with ascending nodes at 120 and 240 degrees, 
respectively . 


The navigation process of GPS proceeds as follows: First, the 

MCS uplinks messages, consisting of time synchronization and SV 
ephemeris information, to the NDSs and the NDSs, in turn, 
continuously broadcast these messages to the user spacecraft. 
Subsequently, the GPSPAC Receiver/Processor Assembly (R/PA) , 
which is the principle GPSPAC subsystem, records and uses the 
information onboard and processes pseudorange and delta- 
pseudorange observations with an Extended Kalman Filter (EKF) to 
calculate an estimate of the user spacecraft's orbit. (This 
information is retained to be analyzed and compared with the 
Landsat definitive ephemeris tape files which are derived from 
independent sources) . If no SVs are in view of the user 

spacecraft, the R/PA of the user spacecraft must propagate its 
own orbit by using a numerical integrator. 


One aspect of the GPSPAC experiment was to validate and to 
calibrate the accuracy of the orbit information produced by GPS 
data; another aspect of the GPSPAC experiment was to determine 
ways to improve the GPSPAC Kalman Filter's navigation performance 
by investigating various data base constant changes or by 
adopting algorithmic changes to the GPSPAC software. To support 
these efforts, a ground-support-modular system called the 
GPSPAC/Landsat-D Interface (GLI) System was developed in March 
1982 by GSFC's Systems Development Branch. The GLI system 
consists of five subsystems and the function of each subsystem is 
described below. 
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First, the GPSPAC Experiment Data Preprocessor (GEDAP) is the 
front-end of the system; it reads, sorts, and reformats telemetry 
strip tape files containing raw GPSPAC measurement data 
(observations and residuals) and GPSPAC Kalman Filter parameter 
estimates. Second, the COMPAR subsystem compares ephemeris files 
from different sources. Next, the PLOT subsystem generates 
graphs of observations, residuals, and filter parameter estimates 
as well as ephemeris differences from COMPARE. Finally, the last 
two subsystems, RECON and ONPAC, are used sequentially. RECON 
recombines files from GEDAP output to produce data in a form 
ready for ONPAC to use. ONPAC, the onboard navigation package, 
is a menu-driven system which has two functions: (1) It uses the 
recombined files generated from RECON to produce estimates of the 
GPSPAC navigation solutions by emulating the GPSPAC Kalman Filter 
data processing scheme. (2) It simulates estimates of the 
GPSPAC navigation solutions by allowing the user to change 
various filter parameters; therefore, one can analyze the effect 
on the GPSPAC navigation solutions once the options are invoked. 


The telemetry strip tape files and the Landsat definitive 
ephemeris tape files, which the GLI system processes, are 
provided by the Landsat Operations Control Center and by the 
Ground Spaceflight Tracking and Data Network (GSTDN) Center, 
respectively. The Landsat Operations Control Center, located at 
GSFC Building 28, retrieves the telemetry strip tape information 
from the playback recordings of the GPSPAC during the satellite 
flyby of a ground tracking station. The GSTDN center, located at 
GSFC Building 25, collects Unified S-Band (USB) range and range- 
rate data. Then, the Goddard Trajectory Determination System 
(GTDS ) is used to process the GSTDN USB data and to compute the 
definitive orbits by performing batch-least squares orbit fits. 


3.0 ONPAC AND GPSPAC ANOMALIES 


From March 1982 until August 1986, the GLI system was operated 
successfully to give the orbit information needed to validate the 
GPSPAC. Namely, comparisons with definitive ephemeris indicated 
that errors in Landsat-4 and 5 position and velocity from GPSPAC 
were consistently less than 50 meters and 6 cm/sec, respectively, 
during periods of good NDS SV visibility (generally speaking 4 
SVs in view) , and that the peak position errors were generally 
less than 1,500 meters during periods of poor SV visibility. 
Although the GLI system helped us to assess the validity of the 
GPSPAC, it has also enabled us to detect two anomalies: one 
pertains to the GLI system and the other pertains to the pre- 
operational phase of GPS. A description of these anomalies is 
given in the next subsections. 
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3 . 1 The ONPAC Anomaly 


A problem which pertains to the ONPAC system is depicted by the 
graphs of the ONPAC position (velocity) uncertainty parameters; 
these graphs are inconsistent with the graphs of the GPSPAC 
position (velocity) uncertainty parameters. For instance, 
figures 1 and 2 illustrate the inconsistencies between the two 
position uncertainty parameters where each graph was generated 
during the same arbitrary time span. These inconsistencies 
suggest that possibly the ONPAC orbit propagator does not match 
the models of GPSPAC entirely. To understand why the ONPAC 
position uncertainty parameters graphs differ from the GPSPAC 
position uncertainty parameters graphs, the ONPAC position 
differences graphs were analyzed; to resolve the ONPAC 
inconsistencies, the ONPAC software was investigated and 
modified. Briefly, the steps taken to resolve the ONPAC 
inconsistencies involved comparing the GPSPAC navigation code 
against the supposedly equivalent ONPAC code, modifying the non- 
conforming routines, and comparing hand-calculated values of 
various filter parameters (based on the GPSPAC algorithms) with 
the values used by ONPAC. A detailed explanation of these 
methods is given in section 4.1. 


3 . 2 The GPSPAC Anomaly 


A problem which is characteristic of the pre-operational phase of 
GPS, is that the GPSPAC/GSTDN definitive position (velocity) 
differences tend to fluctuate tremendously, when the user 
spacecraft is forced to propagate its own orbit because of poor 
SV visibility. For instance, figure 3 illustrates a typical 
graph of the GPSPAC/GSTDN definitive position difference 
fluctuations. Also, Figure 4 shows the NDS visibility to the 
Landsat-4 spacecraft for that period ( only NDS SV numbers 5,6,8 
& 9 were operational for that period) . Notice how the 
GPSPAC/GSTDN definitive position differences graph peaks whenever 
there are less than two SVs in view during any particular time 
span. These fluctuations suggest that there could have been some 
inconsistencies between the way that the GPSPAC orbit propagator 
was designed and implemented. Therefore, to investigate this 
suggestion fully, the GPSPAC navigation software design was 
compared and analyzed with the software code. In addition to 
this, ONPAC was used to simulate runs of the GPSPAC definitive 
position difference fluctuations which enabled us to recommend 
ways to reduce the actual GPSPAC fluctuations. Briefly, the 
steps taken to reduce the simulated GPSPAC fluctuations involved 
studying the GPSPAC navigation design, comparing the design with 
the actual code (to see if the formulas were implemented 
correctly) , documenting the differences, and changing various 
ONPAC filter parameters. A detailed explanation of these methods 
is given in section 4.2. 
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Figure 1. ONPAC Position Uncertainty Parameters 
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Figure 2. GPSPAC Position Uncertainty Parameters 
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NDS VISIBILITY RSS POSITION DIFFERENCES IN METERS 
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Figure 3. GPSPAC/GSTDN Definitive Position Difference Fluctuations 
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Figure 4. NDS Visibility to Landsat-4 
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4.0 STEPS TAKEN TO RESOLVE /REDUCE THE ANOMALIES 

4.1 Steps Taken to Resolve the ONPAC Anomaly 


Since the ONPAC orbit propagator did not match the orbit 
propagation model of GPSPAC entirely, the ONPAC software was 
compared with the GPSPAC software to see where the differences 
occurred. Moreover, several of the filter parameters were 
calculated by hand to check the computed answers given by ONPAC. 
A brief description of the GPSPAC/ONPAC code differences, the 
filter parameters calculated, and the implications of the changes 
is given in the next subsections. 


The GPSPAC/ONPAC Code Differences 


The GPSPAC/ONPAC code differences that were found by comparing 
the two software systems were minor and the following 
modifications were made to the ONPAC code: First, the 

geopotential force model in ONPAC, which is one of the modeled 
external forces used to describe the equations of motion for 
orbit propagation, was upgraded from a 4x4 earth geopotential 
model to a 5x5 earth geopotential model to match that of GPSPAC. 
Next, the atmospheric density model in ONPAC, which is used to 
model the external drag force (another external force used in 
orbit propagation) , was assigned the same lowest altitude 
threshold value as that of GPSPAC. Finally, a variable used in 
ONPAC to validate the pseudorange observations, was replaced with 
another variable to help simulate the measurement error 
computation better. 


The Filter Parameters Calculated 


Another vehicle used to help locate the GPSPAC/ONPAC code 
differences was to calculate by hand the following EKF parameters 
given the GPSPAC EKF software, an arbitrary state vector, and the 
corresponding state-error covariance matrix (see the heading 
entitled "The EKF Background" for a detailed explanation of these 
EKF parameters) : (1) The arbitrary state vector and the state- 
error covariance matrix were propagated to the pseudorange 
measurement time. (2) The pseudorange measurement residual and 
Kalman Gain were calculated. (3) The state vector was updated by 
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Figure 5. ONPAC Position Uncertainty Parameters using the modified code 
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4.2 Steps Taken to Reduce the Simulated 6PSPAC Anomaly 


In order to find out why the GPSPAC fluctuations were occurring, 
the design of the GPSPAC software and the actual software code 
were compared to see if there were any inconsistencies between 
the way the navigation scheme was designed and implemented. 
Basically, this involved studying the EKF algorithms to see how 
the navigation solution was propagated and estimated. In order 
to reduce the simulated GPSPAC fluctuations, various filter 
parameter changes were invoked using ONPAC. What follows in the 
next subsections is a description of the EKF (the source of the 
design-code differences) , the results of filter parameter 
changes, and the implications of these steps taken to reduce the 
simulated GPSPAC fluctuations. 


The GPSPAC Design-Code Differences 


The Extended Kalman Filter Background 


The Extended Kalman Filter (EKF) is the essential element of the 
GPSPAC navigation software. The EKF is an algorithm that 
computes an optimal estimate of the state of a non-linear system, 
given measurements, initial conditions, and statistical 
parameters. The filtering algorithm requires two input 
parameters: an estimate of the state at a previous measurement 
time and an estimate of the state-error covariance matrix at a 
previous measurement time. Given the input parameters, the 
filtering process proceeds as follows: 


(1) The previous filter state, defined as a 9-state vector 

where components 1-3 and 5-7 are the current position 
and velocity of the state and components 4, 8, and 9 

are the user spacecraft's receiver clock time bias, 
receiver clock frequency bias, and satellite drag 
factor, respectively, is propagated to the pseudorange 
measurement time (distance from the NDS satellite to the 
user spacecraft divided by the speed of light 
uncorrected for user clock error) . 

(2) The previous filter state-error covariance matrix, 

defined as a 9x9 matrix where the filter state error is 
given on the main diagonal, is propagated to the 
pseudorange measurement time. 
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(3) The pseudorange measurement gain is calculated to 
determine how much can be "gained" from the 
measurement; then the measurement residual is used to 
determine an estimate of the state error at the 
pseudorange measurement time. 

(4) The delta-pseudorange measurement (the difference 
between two pseudorange measurements) gain is 
calculated and the measurement residual is used to 
update the state error estimate of step 3 . 

(5) The updated state error estimate is used to correct the 
propagated filter state from step 1 and the measurement 
gain calculated is used to update the propagated filter 
state-error covariance matrix from step 2. This results 
in a new filter state and a new state-error covariance 
matrix applicable at the pseudorange measurement time. 


Uyeminami describes the EKF process in detail (2) . 

Figure 6 illustrates the five steps of the filtering process. 



Figure 6. R/PA Extended Kalman Filter 
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The Source of the Design-Code Differences 
The State-Error Covariance Matrix 


The major design-code differences that were found pertained to 
the implemented state-error covariance matrix, P n . The state- 
error covariance matrix at the current measurement time is 
computed as: P n = 7 + Q«D . Briefly, the state 

transition matrix, 31 , propagates the state-error covariance 

matrix forward in time and the state-process-noise covariance 
matrix, Q n , is computed to compensate for the neglected terms 
in the force model of the state. See the Mathematical 
Specifications of the Onboard Navigation Package (ONPAC) 
Simulator for a detailed description of the state-error 
covariance matrix derivation (3) . 


Now, the only way to verify that the software was computing the 
elements of the state-error covariance matrix exactly the way the 
design plans had specified was to check and see if the code was 
computing the elements of the component matrices properly-- 
namely, ^5 and Q_. So, we inspected the equations in the 
software for both tne 3E and Q n matrices and computed, by hand, 
the 81 elements in each matrix; as a result, we discovered that 
the design plan was inconsistent with the code for both matrices. 

The Matrix Design-Code Differences 

Nine elements near the bottom left-hand side of the matrix 

pertain to the modeled geopotential acceleration; these elements 
appeared in the design plans of the matrix, yet, for reasons 
unknown, they were not computed in the GPSPAC software code (see 
figure 7). Consequently, the equations of acceleration due to 
the geopotential were not modeled in the GPSPAC software's 
version of the estimated corrections to the state. 

The Q n Matrix Design-Code Differences 
The design plans for the Q n matrix is the following: 
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Figure 7. State Transition Matrix 
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I = 3x3 Identity Matrix 

Qb = The user spacecraft's process noise time bias term 

Q6 = The user spacecraft's process noise frequency bias term 

Qbb = The user spacecraft's process noise time bias/frequency 

bias coupling term 


This formulation includes several terms that pertain to a modeled 
rotational force? these terms are underlined above. They were 
not computed in the GPSPAC software code because, during the 
design-code phase, the magnitude of these terms were judged to be 
insignificant by the design team. Consequently, the rotational 
force was not modeled in the GPSPAC software's version of the 
state-process-noise covariance matrix. 
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The Results of the Filter Parameter Changes 


Heuberger used Landsat-4 data and ONPAC to analyze the effect on 
various navigation solutions by invoking the option to reduce the 
integration step size, by invoking the option to expand the state 
transition matrix eg and by invoking the option to "tune" the 
state-process-noise covariance matrix Q n . He discovered that by 
incorporating these changes before a simulated run, the simulated 
GPSPAC fluctuations were reduced considerably. Consequently, he 
recommended that these changes should be incorporated into the 
GPSPAC Kalman Filter software. For a detailed explanation of his 
results, see his paper entitled, "The Landsat-4/GPS Experiment 
Final Report". 


A follow-up study was done with Landsat-5 data to test 
Heuberger' s results and to reinforce his recommendations. This 
was accomplished by studying 2 specific arcs from the Landsat-5 
data collection and by using ONPAC to invoke the same filter 
parameter options discussed above. These results, which concur 
with Heuberger' s, are presented below in tabular and graphical 
form; also, a description of the filter parameters invoked is 
provided. 


TABLE 1 

GSTDN DEFINITIVE EPHEMERIS MAXIMUM ERRORS 
OVER SELECTED 10 -HOUR DATA ARCS 


Arc i 

Start Time 

Position (m) 

Velocity (cm/sec) 

1 

June 12, 1984, 05 h 

48.7 

4.2 

2 

July 24, 1985, 03 h 

40.4 

3.3 


In Table 1 the two selected 10-hour data arcs are defined by 
their start dates and times. Furthermore, the maximum 

position/velocity errors of the GSTDN definitive ephemeris are 
given over each of the 10-hour data arcs. The maximum 

position/velocity errors were obtained by performing orbit fits 
over 24 to 32-hour tracking arcs with some overlap between 
successive arcs. The comparison of the two sets of GSTDN 
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definitive ephemeris over the common time span gives one a 
measure of the consistency between the two orbit solutions. It 
was assumed that the maximum error in position (velocity) over 
the two definitive arcs was less than the maximum position 
(velocity) difference in the overlap if the RMS of the USB range 
(range-rate) residuals from the two least squares fit was small 
(i.e. if the data fits were good) . 


TABLE 2 

GPSPAC vs. ONPAC NAVIGATION ACCURACY 


Arc # 

Filter Options 

Eohemeris 

Differences 

■ 

Geopotential 
terms in M 

hhh 

Position (m) 
MAX RMS 

Velocity (cm/sec) 
MAX RMS 

1* 

3 . 0 

No 

CO 

i 

o 

H 

720 

268 

184 

43 

1 

1.0 

No 

CO 

1 

o 

H 

625 

170 

70 

22 

1 

1.0 

Yes 

CO 

1 

o 

H 

526 

165 

65 

22 

1 

1.0 

Yes 

10 6 

466 

160 

65 

21 

2* 

3.0 

No 

00 

1 

o 

H 

539 

128 

83 

19 

2 

1 

No 

10~ 8 

346 

102 

63 

15 

2 

m 

Yes 

CO 

1 

o 

H 

288 

88 

63 

14 

2 

Hi 

Yes 

10 6 

238 

83 

57 

12 


*GPSPAC Solution 


In Table 2, the ONPAC results from invoking the various filter 
parameter changes are summarized. The GSTDN definitive ephemeris 
was compared with the navigation solutions from GPSPAC as well as 
the navigation solutions from ONPAC. The runs were compared 
during the last 6 hours of each data arc to decrease the length 
of time required for ONPAC to process the navigation solutions. 
Moreover, the maximum position/velocity differences and their 
corresponding RMS for each one of the filter parameter changes 
invoked are given for one to analyze and to compare the effect on 
each navigation solution. A plot of GPSPAC/GSTDN definitive 
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position differences over the last 6 hours of arc #1 is shown in 
Figure 8 and the NDS visibility for this period is plotted in 
Figure 9. Similarly, the corresponding plots over the last 6 
hours of arc #2 are shown in Figures 10 and 11. 


By comparing the graphs of the GPSPAC/GSTDN definitive position 
differences, Figures 8 and 10, with the graphs of the NDS 
visibility Figures 9 and 11, one concludes that the GPSPAC 
fluctuations are large due to state error growth during periods 
of prolonged propagation. Also, the fact that GPSPAC uses a 
simple integration scheme— modified Euler with one derivate 
evaluation per step— definitely increases the risk of incurring 
large error growth during poor SV visibility. However, just by 
reducing the step size from 3 seconds to 1 second, the error 
growth became significantly bounded. For the comparison using 
arc #1, refer to Figure 8, the GPSPAC/GSTDN definitive position 
differences; Figure 12, the ONPAC/GSTDN definitive position 
differences with no changes and Figure 13, the ONPAC/GSTDN 
definitive position differences with 1.0 step size. Likewise, 
Figures 10, 16, and 17 give the corresponding comparison for arc 

# 2 . 


A previous section of the paper explains the fact that the design 
plans for the state transition matrix, jiE , included the modeled 
geopotential acceleration terras, yet, for reasons unknown, these 
terms were not included in the code. However, ONPAC offers the 
capability of including the gravity acceleration terms as a 
filter parameter option. So this option was exercised over both 
data arcs and Table 2 shows just how much the maximum 
position/velocity differences decreased by expanding £& to 
include the geopotential acceleration terms. Figures 14 and 18 
illustrate the reduction seen in the GPSPAC fluctuations by using 
a smaller step size, 1.0 second and by using an expanded 3E which 
included the geopotential acceleration terms for each time span. 


A process called "tuning the filter" was exercised over the 2 
arcs in order to compensate for the GPSPAC Kalman Filter's 
underestimate of the true error during periods of poor sv 
visibility. To "tune the filter" one has to adjust the position- 
velocity components of the state-process-noise covariance ^matrix, 
Q n , which are proportional to the date base constant — the 

unmodeled acceleration variance; so by changing Cf£ accordingly, 
the filter is kept from diverging and hence, it produces a better 
estimate of the true error during periods of prolonged 
propagation due to poor SV visibility. The constant was 
increased from 10 -8 m 2 /sec 4 to 10 6 m 2 /sec 4 (see Table 2). This 
filter parameter adjustment improved the error dynamics model 
significantly. Figures 15 and 19 show the ONPAC/GSTDN definitive 
position differences with all three filter parameter options 
invoked; the 1.0 step size reduction, the expanded jjH matrix 
which included the geopotential acceleration terms, and the tuned 
filter result from increasing . 
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Figure 8. GPSPAC/GSTDN Definitive Position Differences - Arc #1 
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Figure 9. NDS Visibility to Landsat-5 - Arc #1 
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Figure 10. GPSPAC/GSTDN Definitive Position Differences - Arc #2 
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Figure 11. NDS Visiblity to Landsat-5 - Arc #2 
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The Implications of the Steps Taken to Reduce the Simulated 

GPSPAC Fluctuations 


The comparison of the GPSPAC navigation design with the code 
proved to be a step in the right direction for two reasons. 
First of all, it helped us detect a couple of problems that 
occurred in the way the EKF state-error covariance matrix was 
designed and implemented— namely, the <I§ and design-code 

differences. Secondly, although the magnitude of the rotational 
force terms of the Q n matrix were indeed insignificant (proven by 
hand calculations) , expanding Jj|j> to include the geopotential 
acceleration terms did help to reduce the simulated GPSPAC 
fluctuations considerably (this result is shown in the section 
entitled, "The Results of the Filter Parameter Changes") . 


In addition to the expanded 3E matrix, two other filter parameter 
options invoked during a ONPAC simulated run also helped to 
reduce the simulated GPSPAC fluctuations: decreasing the 

integration step size and increasing CT<£. to tune the filter. 
Because these results were shown in a follow-up study using 
Landsat-5 data, Heuberger's software recommendations are 
reinforced and are very easy to accommodate — the step size can 
easily be reduced to 1.0 second; the 3? matrix can easily be 
expanded to include the geopotential acceleration terms ; and the 
Q n matrix can easily be tuned by just changing the data base 
constant, CT<£.* 


5.0 CONCLUSIONS 


A GPSPAC/Landsat-D Interface (GLI) ground support system was 
built to validate the performance and to calibrate the accuracy 
of the experimental navigation package, GPSPAC, flown on the 
Landsat-4 and 5 spacecraft. Although the GLI system has operated 
successfully to give the orbit information needed to validate the 
GPSPAC, it also detected the following two anomalies. The first 
problem which pertains to the ONPAC system is that the ONPAC 
orbit propagator is inconsistent with the orbit propagation model 
of GPSPAC. The second problem, which pertains to the pre- 
operational phase of GPS, is that the GPSPAC position (velocity) 
difference fluctuations are large whenever the user spacecraft is 
forced to propagate its own orbit because of poor SV visibility. 
Two attempts were made to resolve the ONPAC inconsistencies: (1) 
comparing the GPSPAC navigation code against the supposedly 
equivalent ONPAC code and modifying the non-conforming routines. 
(2) hand-calculating various filter parameters (to see if these 
answers matched the answers given by ONPAC) . However, neither 
one of them helped; but, it turned out that by invoking the same 
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filter parameter options which were used to reduce the simulated 
GPSPAC fluctuations also helped to reduce the ONPAC 
inconsistencies. The steps taken to reduce the simulated GPSPAC 
fluctuations were discovered by Heuberger. He recommended that 
the following changes should be made to the GPSPAC software: (1) 
reducing the integration step size from 3.0 seconds, to 1.0 
second, (2) expanding the state transition matrix to include the 
geopotential acceleration terms^ (3) increasing the unmodeled 
acceleration variance to tune the Extended Kalman Filter. A 
follow-up study using ONPAC and Landsat-5 data was done to test 
Heuberger's recommended changes. The results of the study 
concurred with his; therefore, his recommendations are 
reinforced. 


Theoretically, several possible software changes could be made to 
ONPAC to reduce the simulated GPSPAC fluctuations, such as 
upgrading the numerical integration scheme to a Runge-Kutta 
integration scheme. But, the software would have to be re- 
designed, re-built, and re-tested which would be costly. 
However, the recommendations mentioned above can be easily 
incorporated in the GPSPAC software thereby making them the 
preferred solution. 


The future of GPSPAC remains to be seen since it is uncertain if 
Landsat or any other spacecraft, for that matter, will ever fly 
another GPS navigation package. However, since August 1986, the 
Systems Development Branch was finished with the project in terms 
of collecting, processing, and analyzing GPSPAC data; we will 
always keep a consolidated collection of good GPSPAC continuous 
data arcs from the Landsat-4 and 5 spacecraft for future 
independent studies of autonomous onboard navigation systems — for 
example, there is speculation that GPS may be used on Space 
Station and we will remain as a point of contact for obtaining 
various information pertaining to the GPSPAC/Landsat-D Interface 
ground support system. 
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